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A SIMPLE MODEL FOR THE ABSORPTION OF STARLIGHT BY DUST IN GALAXIES 
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ABSTRACT 

We present a new model to compute the effects of dust on the integrated spectral properties of galaxies, 
based on an idealized prescription of the main features of the interstellar medium (ISM). The model includes 
the ionization of Hn regions in the interiors of the dense clouds in which stars form and the influence of the 
finite lifetime of these clouds on the absorption of radiation. We compute the production of emission lines and 
the absorption of continuum radiation in the Hn regions and the subsequent transfer of line and continuum 
radiation in the surrounding Hi regions and the ambient ISM. This enables us to interpret simultaneously all 
the observations of a homogeneous sample of nearby ultraviolet-selected starburst galaxies, including the ratio 
of far-infrared to ultraviolet luminosities, the ratio of Ha to H/3 luminosities, the Ha equivalent width, and the 
ultraviolet spectral slope. We show that the finite lifetime of stellar birth clouds is a key ingredient to resolve an 
apparent discrepancy between the attenuation of line and continuum photons in starburst galaxies. In addition, 
we find that an effective absorption curve proportional to A -0 ' 7 reproduces the observed relation between the 
ratio of far-infrared to ultraviolet luminosities and the ultraviolet spectral slope. We interpret this relation most 
simply as a sequence in the overall dust content of the galaxies. The shallow wavelength dependence of the 
effective absorption curve is compatible with the steepness of known extinction curves if the dust has a patchy 
distribution. In particular, we find that a random distribution of discrete clouds with optical depths similar to 
those in the Milky Way provides a consistent interpretation of all the observations. A noteworthy outcome of 
our detailed analysis is that the observed mean relations for starburst galaxies can be closely approximated by 
the following simple recipe: use an effective absorption curve proportional to A -0,7 to attenuate the line and 
continuum radiation from each stellar generation, and lower the normalization of the curve typically by a factor 
of 3 after 10 7 yr to account for the dispersal of the birth clouds. This recipe or our full model for absorption can 
be incorporated easily into any population synthesis model. 

Subject headings: galaxies: ISM — galaxies: starburst 
ultraviolet: galaxies 



infrared: galaxies — ISM: dust, extinction 



1. INTRODUCTION 

To interpret the observed spectral properties of galax- 
ies, we require models for both the production of stellar 
radiation and its transfer through the interstellar medium 
(ISM). Currently, the accuracy of population synthesis 
models contrasts with the rudimentary way in which the 
absorption of starlight by dust is often treated. In many 
applications, dust is cither ignored or assumed to be dis- 
tributed in a uniform screen in front of the stars. The 
resulting uncertainties in the absorption of the ultravio- 
let radiation in galaxies can be as much as an order of 
magnitude. This problem has become especially acute 
in studies of galaxies at high redshifts. Nearby starburst 
galaxies may be suitable analogs of high-redshift galaxies 
and provide important clues to interpreting their spectral 
properties. The observations of nearby ultraviolet-selected 
starburst galaxies are numerous, including the ratio of far- 
infrared to ultraviolet luminosities, the ratio of Ha to H/3 
luminosities, the Ha equivalent width, and the ultraviolet 
spectral slope. In particular, there is a remarkably tight 
correlation between far-infrared luminosity and ultravio- 
let spectral slope (Meurer et al. 1995; Meurer, Heckman, 
& Calzetti 1999). This wealth of observations can po- 



tentially help us quantify the effects of dust on various 
spectral properties of galaxies. 

There have been several analyses of the spatial distri- 
bution and optical properties of the dust in nearby star- 
burst galaxies based on various subsets of the observations 
(e.g., Fanelli, O'Connell, & Thuan 1988; Calzetti, Kinney, 
& Storchi-Bergmann 1994, 1996; Puxley k, Brand 1994; 
Meurer et al. 1995; Gordon, Calzetti, & Witt 1997). A 
generic result of these studies is that if the dust is dis- 
tributed in a uniform foreground screen, it must have an 
unusually grey extinction curve. Otherwise, the distribu- 
tion must be patchy. However, these analyses also raise 
questions of how to account self-consistently for all the 
observations. For example, the absorption inferred from 
the Ha/H/3 ratio in starburst galaxies is typically twice as 
high as that inferred from the ultraviolet spectral slope. 
One interpretation of this result is that the ionized gas and 
ultraviolet-bright stars have different spatial distributions 
(Calzetti 1997). Issues such as this highlight the need for 
a simple yet versatile model to interpret simultaneously 
a wide range of phenomena related to the absorption of 
starlight by dust in galaxies. The purpose of this paper is 
to present such a model. 
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We begin with the conventional view that young stars 
ionize Hn regions in the interiors of the dense clouds in 
which they are born. Line photons produced in the Hn re- 
gions and the non-ionizing continuum photons from young 
stars are absorbed in the same way by dust in the outer Hi 
envelopes of the birth clouds and the ambient ISM. The 
birth clouds, however, have finite lifetimes. Thus, non- 
ionizing ultraviolet and optical photons from stars that 
live longer than the birth clouds are absorbed only by 
the ambient ISM. This allows the ultraviolet continuum 
to be less attenuated than the emission lines. Our model 
builds on several previous studies. For example, Silva et 
al. (1998) considered the effects of finite lifetimes of stellar 
birth clouds on the continuum but not the line emission 
from galaxies. Here, we treat the transfer of radiation (es- 
pecially scattering) in an approximate way, which would 
preclude a detailed description of surface brightnesses but 
should be appropriate for angle-averaged quantities such 
as luminosities. Our model succeeds in accounting quanti- 
tatively for all the available observations of a homogeneous 
sample of nearby starburst galaxies. 

We present our model in §2, where we express the effec- 
tive absorption curve describing the global transmission of 
radiation in terms of the different components of the ISM. 
In §3, we compare our model with observations and iden- 
tify the specific influence of each parameter on the different 
integrated spectral properties of galaxies. One outcome 
of our detailed analysis is a remarkably simple recipe for 
absorption, which provides a good approximation to the 
observed mean relations and is easy to incorporate into 
any population synthesis model. In §4, we explore how 
the spatial distribution of dust can be constrained by the 
observations. Our conclusions are summarized in §5. 

2. DEFINITION OF THE MODEL 

We first introduce some notation and nomenclature that 
will help specify our model. The luminosity per unit wave- 
length L\(t) emerging at the time t from a galaxy illumi- 
nated by an internal stellar population can be expressed 
generally as 



by the luminosity of each stellar generation T\(t) and the 
total unattenuated stellar luminosity: 



L x (t) = T x (t) / dt'^{t-t')S x (t'). 



(2) 



_ ^ J*dt>iP(t-t>)S x (t>)T x (t,t>) 
f*dt>iP(t-t>)S x (t>) 

We often reexpress T x (t) in terms of an "effective absorp- 
tion" optical depth f A (i) indicated by a caret and defined 
by 

T A (i)=exp[-f A (t)]. (4) 

This nomenclature emphasizes that the luminosity of a 
galaxy is diminished purely by the absorption of photons. 2 
However, since scattering may also affect the path lengths 
of the photons through the ISM before they are absorbed, 
r x (t) will differ in general from the true absorption curve 
determined by the optical properties of the dust grains. 
The fraction of stellar radiation not transmitted by the 
ISM, 1 — T x (t), is absorbed by dust. Using equation (2) 
and integrating over all wavelengths, we obtain the total 
luminosity absorbed and reradiated by dust 

POO pt 

L dust (t)= d\[l-T x (t)] / dt'iP(t-t')S x (t'). 
Jo Jo 



Lx(t) 



dt'iP(t-t')S x (t')T x (t,t'). 



(1) 



( 5 ) 

We now specify our model for computing the above 
properties of galaxies. This amounts to specifying the 
transmission function T x (t,t r ). Our goal is to parame- 
terize only those features of the ISM that are essential 
to understanding the main physical effects of dust on the 
integrated spectral properties of galaxies. We adopt the 
simple but fairly conventional view illustrated in Figure 1. 
Young stars ionize Hn regions in the inner parts of the 
dense (mostly molecular) clouds in which they are born. 
Line photons produced in the Hn regions and the non- 
ionizing continuum photons from young stars propagate 
through the outer Hi envelopes of these "birth clouds" 
and then through the "ambient ISM" before they escape 
from the galaxy. Both the Hn and Hi regions of the birth 
clouds may contain dust, and in the ambient ISM, the 
dust may be distributed in a smooth or patchy way. The 
birth clouds, however, have finite lifetimes. Thus, photons 
emitted by longer-lived stars propagate only through the 
ambient ISM — where they are treated in the same way 
as photons emerging from the birth clouds — before they 
escape from the galaxy. To incorporate these basic fea- 
tures into our model, we make the following simplifying 
assumptions. 

We assume that stars form in clusters with the galaxy- 
wide stellar initial mass function (IMF). This should be 
an excellent approximation for the purpose of computing 
the integrated spectral properties of a galaxy, although it 
could break down in clusters or regions with small num- 
bers of stars. In our model, all star clusters are embedded 
in their birth clouds for some time and then disrupt them 
or migrate away from them into the ambient ISM of the 
galaxy. We assume for simplicity that the effective absorp- 
tion in the birth clouds depends only on the stellar age t' 

2 The optical depth r x defined by equation (4) has also been referred to by several other names in the literature, including the optical depth 
of "apparent extinction," "effective extinction," "attenuation," "obscuration," and "absorption." 



Here ip(t — t') is the star formation rate at the time t — t', 
S\(t') is the luminosity emitted per unit wavelength and 
per unit mass by a stellar generation of age t', and T x (t, t') 
is the "transmission function," defined as the fraction of 
the radiation produced at the wavelength A at the time t by 
a generation of stars of age t' that escapes from the galaxy. 
Thus, T x (t,t') must be regarded as the average transmis- 
sion along rays emanating in all directions from all stars of 
age t' within the galaxy. Since absorption and scattering 
by dust will cause the radiation to emerge anisotropically 
to some degree, L x can only be measured by a hypotheti- 
cal set of detectors completely surrounding the galaxy. In 
practice, this means that equation (1) gives the mean of 
AttD 2 /a, where D is the distance and f x is the observed 
flux, for an ensemble of randomly oriented but otherwise 
similar galaxies. 

For some purposes, it is convenient to reexpress L x (t) 
as the product of a mean transmission function weighted 
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and is constant in the ambient ISM. Hence, we can write 
the total transmission function in the form 



T x (t,t') = T* c (t')T l x 



(6) 



where c and Tj[ SM are the transmission functions of the 
birth clouds and the ambient ISM, respectively. In our 
model, the birth clouds, assumed to be identical, consist 
of an inner Hll region ionized by young stars and bounded 
by an outer Hi region. For simplicity, we assume that 
all stars remain within their birth clouds with a constant 
transmission function for a time iec and are then intro- 
duced immediately into the ambient ISM; thus 



7f c (0 = 



for t' < t B c , 
for t' > t BC , 



(7) 



where T™ and T™ are the transmission functions of the 
Hn and Hi regions, respectively. By writing the transmis- 
sion functions as products in equations (6) and (7), we 
have assumed that the propagation of photons in the Hn 
and Hi regions and the ambient ISM are independent of 
each other. This prescription is exact for forward scatter- 
ing and should be a good approximation for the purpose 
of computing luminosities even with some backward scat- 
tering. 

An alternative but equivalent way of expressing the 
transmission functions of the Hn and Hi regions and the 
ambient ISM is in terms of the corresponding effective ab- 
sorption optical depths, t™, t™, and f A SM . For non- 
ionizing photons (with A > Al = 912 A), we thus write 



rriHII 
1 X 


= exp(-f A Hn ), 


(8) 


rpHI 

1 A 


= cxp(-f A HI ) , 


(9) 


^ISM 
J A 


= exp(-ff M ) . 


(10) 



For some purposes, it is also convenient to define the to- 
tal effective absorption optical depth of dust in the birth 
clouds and the fraction of this contributed by the internal 
Hn regions 

r A BC = r™ + i!", (11) 



We note that / may or may not coincide with the ratio of 
the Hn to total gas column densities of the birth clouds, 
depending on whether or not the dust-to-gas ratios in the 
Hn regions are the same as those in the Hi regions. 

Finally, we assume that the emission lines are produced 
only in the Hll regions of the birth clouds. The justifica- 
tion for this is that the lifetimes of the birth clouds are 
generally greater than the lifetimes of the stars producing 
most of the ionizing photons, i.e. about 3 x 10 6 yr. In this 
case, line photons are absorbed in the Hi regions and the 
ambient ISM in the same way as non-ionizing continuum 
photons. The emergent luminosity of a line of wavelength 
A; is thus 

Li(t) = dt' W - t') Stit'^Tl™ > (13) 

Jo 

where Si(t') is the luminosity of the line produced in the 
Hn regions by a stellar generation of age t' , and T™ and 
yiSM are gj ven by equations (9) and (10). The luminosities 
of H-recombination lines are proportional to the ionization 
rate, 

jfL dxxs " (14) 



N- — 

1 v ion — 



where T- mn is the fraction of Lyman continuum photons ab- 
sorbed by the gas rather than the dust in the Hn regions. 
We adopt the formula 



( fH")3exp(-fH") 



Ti ° n 3{(fH")2 - 2fH" + 2[1 - exp(-fH")]} 



(15) 



(12) 



This equation, which is equivalent to equation (5-29) and 
Table 5.4 of Spitzcr (1978), was derived by Petrosian, Silk, 
& Field (1972) and shown by them to provide a good ap- 
proximation to the results of more detailed numerical com- 
putations. 

We are especially interested in the emergent luminosities 
of the Ha and H/3 Balmer lines of hydrogen, L Ba and Lh/3- 
We assume that, given A^ on , the production rates of all H- 
recombination photons are the same as for the dust-free 
case B recombination. This is validated by detailed pho- 
toionization models of Hn regions, including the effects of 
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0.450 {he/ Ana) N ion , 
0.117 (hc/\ H p)N ion , 



(16) 



internal dust (Hummer & Storey 1992; Ferland 1996; Bot- 
torff et al. 1998). For electronic densities n e <C 10 4 cm~ 3 
and temperatures T e w 10 4 K, we then have (Osterbrock 
1989) 

Sup 

with Ahq = 6563 A and Ah/j = 4861 A. It is straight- 
forward to extend this procedure to compute the lumi- 
nosity of any other H-recombination line (e.g., P/3, B7). 
We assume that all the power radiated in the form of 
La photons, S^a = 0.676 (hc/X^a) A; on , is eventually ab- 
sorbed by dust as a consequence of resonant scattering 
(Alo = 1216 A). We have checked that our conclusions are 
not affected even if a significant fraction of La photons 
escapes from the galaxy. For completeness, we compute 
the power from ionizing radiation that is neither absorbed 
by dust nor emerges in the Ha and H/3 lines, 

rX L 



Sx=[T i 



dX S\ 



Sha + Su a + Sftfj ) , (17) 



and distribute this uniformly in wavelength between 3000 
and 6000 A (a range that includes most of the relevant 
emission lines). This is only a few percent of all the stellar 
radiation and has a negligible influence on our results. 

We now consider how the transmission function of the 
ambient ISM depends on the spatial distribution and opti- 
cal properties of the dust. A particularly enlightening way 
to express T\ SM is in terms of the probability density for 
the absorption of photons emitted in all directions by all 
stars within a galaxy. We define t\ as the optical depth of 
absorption along the path of a photon, from its emission 
by a star to its escape from the galaxy, including the influ- 
ence of scattering on the path. From now on, we refer to 
t\ simply as the "optical depth" of the dust. 3 We further 
define p(T\)d,T\ as the probability that the optical depth 
lies between t\ and t\ + dr\. Integrating over rays ema- 
nating in all directions from all stars within the galaxy, we 
can then write the transmission function of the ambient 
ISM as 



rpISM 

1 A 



dTAp(TA)exp(-T A ) 



(18) 



If the spatial distribution of the dust is specified, this de- 
termines p(t\), and T\ SM then follows from equation (18). 
For some purposes, however, it is more informative to spec- 
ify p(t\) itself rather than any of the spatial distributions 
of the dust that could give rise to it when evaluating Tj[ SM . 

The optical depth t\ of the dust depends on a combina- 
tion of absorption and scattering. If scattering were only 
in the forward direction, t\ would equal the true absorp- 
tion optical depth. In the case of isotropic scattering, t\ is 
given by the geometrical mean of the true absorption and 
extinction (absorption plus scattering) optical depths, i.e., 

rx = [t~\ (tx + rl)] 1/2 = T a x (l u x )- 1/2 , (19) 

where uux is the single-scattering albedo (Rybicki & Light- 
man 1979). 4 Although this formula was derived for an in- 
finite homogeneous medium, it also provides an excellent 
approximation for more complicated geometries (see, e.g., 



Fig. 5 of Silva et al. 1998, which compares this approxima- 
tion with a Monte Carlo calculation by Witt, Thronson, 
& Capuano 1992). In reality, both observations and mod- 
els suggest that interstellar dust in the Milky Way has 
pronounced forward scattering at ultraviolet and optical 
wavelengths (e.g., Drainc & Lee 1984; Kim, Martin, & 
Hendry 1994; Calzetti et al. 1995). Thus, we expect the 
truth to lie between the idealized formulae tx = and 
tx = t£(1 — ujx)^ 1 ^ 2 - In our model, the details of scat- 
tering and absorption are subsumed in our approximation 
of the optical depth of the dust as a power law in wave- 
length, rx oc A~ m (§4). This approach would be question- 
able if we wished to compute directional quantities such 
as surface brightnesses. However, since we compute only 
angle-averaged quantities (i.e., luminosities), our simple 
treatment of scattering should be an adequate approxima- 
tion. The reason for this is that, since many of the photons 
scattered out of a line of sight will be replaced by other 
photons scattered into it, inaccuracies in the treatment of 
scattering will partially (but not totally) cancel out when 
the emission is averaged over all directions. 

We now consider three specific, idealized distributions 
of dust in the ambient ISM: a uniform foreground screen, 
a mixed slab, and discrete clouds. These are fairly stan- 
dard representations, which may be regarded as illustra- 
tive of a wide range of possibilities. Other models could be 
constructed for more realistic and hence more complicated 
distributions of dust, for example by multiplying the trans- 
mission functions derived from different p(t\)'s. It should 
be recalled in any case that the total transmission function 
in our model is the product of that for the birth clouds and 
that for the ambient ISM, as specified by equation (6). 

Foreground screen . The simplest model is that of a uni- 
form foreground screen, which absorbs the light from all 
stars by the same amount. While unrealistic, this model 
has been adopted in many previous studies, and we include 
it here for reference. We call t^ c the optical depth of the 
screen. We then have 



nISM 



S{tx-t?), 
-- exp(-r! c ) , 



(20) 
(21) 



where 5 is the Dirac delta function. 

Mixed slab. A potentially more realistic model is that 
of a slab filled with a uniform mixture of sources and ab- 
sorbers (i.e., stars and dust). This can account for differ- 
ences in the absorption affecting stars near the midplanc 
or near the surface of a galactic disk. We call r|' the full 
optical depth through the slab in the direction normal to 
the surface. By considering sources located at all points 
within the slab and rays emanating at all angles from the 
normal, we obtain, after some manipulation, 



, , fl/(2r| 1 ), 



for tx < T^ 
for tx > r? 1 



(22) 



J.ISM _ 



= [1 + (rf 1) cxp(-rl 1 ) - {iffE^-rf)] 



where E\ is the exponential integral of the first order. 



(23) 



3 The optical depth denoted here by t\ is sometimes referred to as an effective optical depth and denoted by r» (Rybicki & Lightman 1979). 
In this paper, we use the term "effective optical depth" to refer to the geometry-dependent optical depth t\ defined by equation (4). 

4 Equation (19) is exact in the limit of large extinction optical depths but overestimates t\ in the opposite limit. 
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Discrete clouds . Both the foreground screen and mixed 
slab models are based on smooth representations of the 
ambient ISM. In real galaxies, however, the ISM is orga- 
nized into clouds and other structures. The photons emit- 
ted in different directions by different stars can encounter 
any number of these clouds on their way out of the galaxy. 
Thus, p{t\) may be better approximated by a Poisson dis- 
tribution (Spitzer 1978; Natta & Panagia 1984). We call 
t x the optical depth per cloud, and n the mean number 
of clouds encountered along different lines of sight. The 
mean optical depth is thus t\ — nr^. In this case, we 
have 5 



n=0 



exp(-n)<5(r A - nrj) , 



= exp 



-n(l 



(24) 



(25) 



In the limits r? 



and T\ 



oo, 



we obtain the famil- 
iar expressions Tj[ SM — > cxp(— t\) and T^ SM — ► exp(— n). 
This model should provide an adequate description of the 
absorption along individual random lines of sight, as re- 
quired by equation (18), but it neglects possible correla- 
tions between the absorption along multiple nearby lines 
of sight, which may be important in other applications. 

In principle, we could compute models with arbitrary 
star formation rates. For simplicity, in all the following 
applications, we take the star formation rate to be con- 
stant and include the age t of the stellar population as a 
parameter. This should be regarded as the effective age of 
the most recent burst of star formation, which can range 
from a few times 10 7 yr for a galaxy in an active starburst 
to several times 10 9 yr for a more quiescent galaxy. With 
the assumption ip = const, equations (3), (4), and (6)— (11) 
allow us to reexpress the effective absorption optical depth 
at the time t in terms of the effective absorption optical 
depths in the birth clouds and the ISM as 

(26) 



r x (t)=f^+r^ + Ar x (t), 



with 
Af A (i) 

0, 



-In 



/ tBC dt' S x (t') + e^ 



J t ' BC dt' S\(t') 



/„<**' S x (t>) 



for t <t B c, 
for t>t B c- 
(27) 

We compute S\(t') using the latest version of the Bruzual 
& Chariot (1993) population synthesis code. Our conclu- 
sions would not be affected by the use of other standard 
codes. Unless otherwise indicated, we adopt solar metal- 
licity and a Salpeter IMF truncated at 0.1 and 100 M Q . 

We now have a completely specified model for the ab- 
sorption of starlight by dust in galaxies. The main ad- 
justable parameters describing the birth clouds are: the 
lifetime of the clouds, t B c', the wavelength dependence and 



normalization of t x c : 



and the fraction of f A c contributed 



model (t a and n). The other main parameter is the effec- 
tive age of the starburst, t. We compute the continuum 
luminosity per unit wavelength emerging from a galaxy us- 
ing equations (1) and (6)-(10). The emergent Ha and H/3 
line luminosities, Ln a (t) and Lnp(i), are obtained from 
equations (9), (10), and (13)— (16) . The ratio of the Ha 
line and continuum luminosities is then the emission equiv- 
alent width of the emergent Ha line, W Ba (t) . Finally, we 
compute the far-infrared luminosity of the galaxy, Ldust(i)> 
using equations (3) and (5). 

In the remainder of this paper, we take two approaches 
to compare our model with observations. As equa- 
tions (10) and (18) indicate, the optical properties and 
spatial distribution of the dust affect observable quantities 
only through the effective absorption curve r A SM . Differ- 
ent combinations of t\ and p(t\), however, could lead to 
the same wavelength dependence of f A SM . Thus, it is of 
primary interest to study how the wavelength dependence 
of f A is constrained by observations. This is our goal in 
§3 below. Then, in §4, we explore whether there are exam- 
ples of optical properties and spatial distributions of the 
dust that satisfy these constraints, recognizing that there 
is no unique solution. 

3. COMPARISONS WITH OBSERVATIONS: THE EFFECTIVE 
ABSORPTION CURVE 

In this section, we use observations of nearby starburst 
galaxies to constrain the wavelength dependence of the ef- 
fective absorption curve f A defined by equation (4) without 
specifying the spatial distribution and optical properties 
of the dust. We also identify the specific influence of each 
parameter in the model on different integrated spectral 
properties of galaxies and investigate the physical origin of 
the observed relationships between these properties. We 
assume for simplicity that the effective absorption curves 
of the birth clouds and the ambient ISM have the same 
power-law form, 

r£ c (A/5500 A)~", 



-ISM 



— Ty SM (A/5500 A)" 



(28) 
(29) 
-£ c and 



by dust in Hn regions, /. For the ambient ISM, the pa- 
rameters depend on whether we adopt a foreground screen 
model (t| c ), a mixed slab model (t a ), or a discrete cloud 

5 Strictly, equation (24) corresponds to a Poisson distribution of identical, face-on screens. For spherical clouds, the factor S(t\ — nr|) 
would have to be replaced by a function with a finite width to account for the different optical depths seen by photons with different impact 
parameters. We ignore this subtlety here. 



where n is a parameter to be determined, and 
fy SM are the normalizing coefficients at 5500 A. 

To calibrate our model, we appeal to a homogeneous 
set of observations of nearby ultraviolet-selected starburst 
galaxies, for which the ultraviolet, far-infrared, Ha, and 
H/3 fluxes, Ha equivalent widths, and ultraviolet spectral 
slopes are available. The sample, compiled by Meurer et 
al. (1999), includes 57 non-Seyfert galaxies from the In- 
ternational Ultraviolet Explorer (IUE) atlas of Kinney et 
al. (1993) classified into one of the starburst categories 
(i.e., starburst nucleus, starburst ring, blue compact dwarf 
galaxy, or blue compact galaxy). This sample has been fur- 
ther restricted by Meurer et al. (1999) to galaxies with op- 
tical angular diameters less than 4', for which most of the 
(usually concentrated) ultraviolet emission could be ob- 
served within the IUE aperture. The galaxies span a wide 
range of morphological types, from Sab to Irr, and a wide 
range of absolute i?-band magnitudes, — 16 < Mb < —22 
(we adopt a Hubble constant H n = 70 kms _1 Mpc _1 ). 
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From the IUE spectra of all galaxies in the sample, Meurer 
et al. (1999) compute the restframe ultraviolet flux at 
1600 A, Fi600j and the slope (3 of the ultraviolet contin- 
uum. The quantity -Fieoo is defined as A/a, where /i6oo is 
the mean flux density over a square passband with a cen- 
tral wavelength of 1600 A and a width of 350 A. The slope 
(3 is determined from a power-law fit of the form f\ oc X 13 to 
the observed ultraviolet spectrum in ten continuum bands 
that avoid strong stellar and interstellar absorption fea- 
tures (in particular, the 2175 A dust feature) in the range 
1268 < A < 2580 A, as defined by Calzetti et al. (1994). 
These ultraviolet properties were derived for galaxies in 
the sample after removing foreground extinction by the 
Milky Way. 

Observations with the Infrared Astronomical Satellite 
(IRAS) are available for 47 galaxies in the sample. To es- 
timate the total far-infrared flux fdust from the observed 
IR A S flux densities, /„(60/im) and /^(lOO/im), we follow 
the prescription of Meurer et al. (1999). They first cal- 
culate the quantity F Fm = 1.26 x 10" 11 [2.58/^(60 /urn) + 
/„(100/im)] ergscm _2 s _1 defined by Helou et al. (1988) 
and then compute Fd us t (-Pbol in their notation) from Ffir 
and /i / (60/Ltm)//i / (100/Ltm) using the bolometric correc- 
tion appropriate for dust with a single temperature and 
an emissivity proportional to frequency v. This prescrip- 
tion agrees reasonably well with observations from near- 
infrared to submillimctcr wavelengths for a small sample 
of galaxies (see Fig. 2 of Meurer et al. 1999). The correc- 
tions from Ffir to Fd us t are relatively small, with a mean 
(-Fdust /-Pfir) ~ 1-4 and a standard deviation of 0.1. Fi- 
nally, for 31 galaxies, the ratio of Ha to H/3 fluxes and the 
Ha equivalent width are available from Storchi-Bcrgmann, 
Kinney, & Challis (1995) and McQuade, Calzetti, & Kin- 
ney (1995), with corrections for stellar absorption by 
Calzetti et al. (1994). Neglecting possible anisotropies, we 
equate flux ratios at all wavelengths to the corresponding 
luminosity ratios. Figures 2, 3, and 4 show Ldust/^1600, 
Lna/Lnp, and Wn a , respectively, as a function of (3 for 
this sample. The typical scatter at fixed f3 amounts to 
a factor of about 2 in £dust/^i600 and Lua/Ln^, but 
nearly a factor of 10 in Wn a - Thus, both £dust/£i600 
and Ln a /Lnp, to a first approximation, are functions of 
(3, while there is little or no correlation between Wu a and 
(3. We note, however, that the scatter even in the first two 
relations appears to be real because it is larger than the 
typical measurement errors (which are shown in the upper 
left panels of Figs. 2-4). 

We now use these observations to constrain the param- 
eters in our model. Our goal is to identify models that ac- 
count for the typical properties, scatter, and trends seen in 
the sample, rather than to interpret in detail the emission 
from individual galaxies. Such an attempt might be of 
limited validity because the predictions of our model, un- 
like the observations, are averaged over angles. 6 We have 
checked, however, that the properties of every galaxy in 
the sample could be reproduced with at least one combi- 
nation of parameters. The influence of each parameter on 
observable quantities can best be explored by keeping all 
other parameters fixed at "standard" values. After some 



experimentation, we adopted: 

t = 3 x 10 8 yr 
i BC = lx 10 7 yr 
f^ C = 1.0 
4 SM = 0.5 
/ = 0.1 
n = 0.7. 



(30) 

While these values are not the results of a rigorous op- 
timization procedure, they do enable the standard model 
to match roughly the observed typical (i.e., median) prop- 
erties of the starburst sample. We note that fy C w 1 is 
typical of the observed optical depths to O-type and su- 
pergiant stars in the Milky Way (Humphreys 1978). The 
dense cores of molecular clouds have higher optical depths, 
but even short-lived stars spend only the first 10%— 20% 
of their lifetimes in these regions (e.g., Mezger & Smith 
1977). Also, tec ~ 10 7 yr is typical of the timescale for 
the dispersal of giant molecular clouds in the Milky Way 
(Blitz & Shu 1980). 

We now compare the predictions of our model with the 
observations described above. Each panel in Figures 2-4 
shows the effect of increasing and decreasing one parame- 
ter with respect to its standard value with the others held 
fixed. We can summarize the role of each parameter as 
follows. 

Effective starburst age . The effect of increasing t is most 
noticeable on the Ha equivalent width, which drops from 
almost 300 A at t = 3 x 10 7 yr to less than 100 A at 
t = 3 x 10 8 yr because of the rising contribution by old 
stars to the continuum near the line (Fig. 4a). Increasing 
t also makes the ultraviolet spectral slope larger because 
old stars have intrinsically red spectra. However, this does 
not affect the Ha/H/3 ratio, since the emission lines are 
produced mainly by stars younger than about 3 x 10 6 yr 
(Fig. 3a). 

Lifetime of the birth clouds . Increasing tec makes both 
(3 and the ratio of far-infrared to ultraviolet luminosities 
larger because the radiation from young blue stars is heav- 
ily absorbed for a longer time (Fig. 26). Changes in tec 
do not affect Ln a / Lh/3 because we have assumed that the 
birth clouds last longer than 3 x 10 6 yr, the lifetime of 
stars producing most of the ionizing photons (§2). How- 
ever, since more of the continuum luminosity near Ha is 
absorbed for larger tec, the Ha equivalent width also in- 
creases (Fig. 46). 

Effective absorption optical depth in the birth clouds. 
Increasing fy C above unity has the remarkable property of 
making Lh«/-^h/3 larger at almost fixed (3 and £dust/-^i600 
(Figs. 2c and 3c). The reason for this is that, while 
Lua/ Luf3 traces the absorption in the Hi envelopes of the 
birth clouds, the ultraviolet radiation from embedded stars 
is so attenuated for fy C > 1 that Ldust saturates, and (3 
and iigoo are dominated by stars in the ambient ISM. 
This has fundamental observational implications (see be- 
low). Figure 4 c also shows that the Ha equivalent width 
is very sensitive to changes in fy C , which affect Lu a . 

Effective absorption optical depth in the ambient ISM. 
Since fy SM controls the absorption of stars that have 
emerged from their birth clouds, it has a major influence 

6 Although some differences are expected between the luminosity ratios and flux ratios due to anisotropic emission from galaxies, these must be 
relatively small or else the relations shown in Figures 2 and 3 would be washed out. The obvious outlier at j3 = —0.76 and log(Ld uat /Li60o) = 2.2 
in Figure 2 could be a galaxy for which some of the ultraviolet emission was missed with IUE (see the discussion in §3.1 in Meurer et al. 1999). 
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Fig. 2. — Ratio of far-infrared to ultraviolet luminosities plotted against ultraviolet spectral slope. The data points are from 
the Meurer et al. (1999) sample discussed in §3 and are repeated in all panels (typical measurement errors are indicated in the 
upper left panel). In each case, the line shows the effect of varying one parameter from the lower end of the range [square) to the 
standard value [circle) to the upper end of the range {triangle), with all other parameters fixed at their standard values (eq. [30]): 
(a) effective starburst ages, t = 3 x 10 7 , 3 x 10 8 , and 3 x 10 9 yr; (b) lifetimes of the birth clouds, tsc = 3 x 10 6 , 1 x 10 7 , and 
3 x 10 7 yr; (c) effective optical depths in the birth clouds, Ty C — 0.0, 1.0, and 2.0; (d) effective optical depths in the ambient ISM, 
^ism _ q o 5^ an( j i q. (g) f rac tions of dust in Hn regions in the birth clouds, / = 0.0, 0.1, and 1.0; and (f) exponents in the 
relation ff M oc A" n , n = 0.4, 0.7, and 1.0. 
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Fig. 3. — Hq/H/3 ratio plotted against ultraviolet spectral slope. The data points are from the Meurer et al. (1999) sample 
discussed in §3 and are repeated in all panels (typical measurement errors are indicated in the upper left panel). In each case, 
the line shows the effect of varying one parameter from the lower end of the range (square) to the standard value (circle) to the 
upper end of the range (triangle), with all other parameters fixed at their standard values (eq. [30]). The models are the same as 
in Fig. 2. 
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on the ultraviolet spectra and hence on /3 and £d U st/£i600 
(Fig. 2d). Chan ges in r^p M also affect the Ha/H/3 ratio, 
although with a strong dependence on [3 (Fig. 3c?). As 
expected, fy M has no effect on Wu a - 

Fraction of dust in the HII regions . Increasing / causes 
more ionizing photons to be absorbed by dust and hence 
reduces the Ha equivalent width considerably (Fig. 4e). 
Dust in Hn regions has a negligible effect on the Ha/H/3 
ratio. However, since we have kept Ty C fixed, increasing 
/ reduces the absorption in the Hi envelopes of the birth 
clouds and hence the Ha/H/3 ratio in Figure 3e. Changes 
in / have no effect on [3 and Li6oo an d a negligible effect 
on idust- 

Exponent of the effective absorption curve . Increasing 
n steepens the effective absorption curve. This increases 
the ratio of far-infrared to ultraviolet luminosities, the 
Ha/H/3 ratio, and the ultraviolet spectral slope and re- 
duces slightly the Ha equivalent width (Figs. 2/-4/). We 
elaborate below on the reason for adopting n = 0.7 as a 
standard value. 

Other parameters . We have also computed models with 
different IMFs and metallicities (not shown). Changes in 
these parameters have no influence on Ln a /Lnp and affect 
[3 and Ldust/-^i600 at the level of less than 20%. The influ- 
ence on the Ha equivalent width is slightly stronger. For 
example, adopting a Scalo IMF instead of a Salpeter IMF 
would reduce Wu a by about 50% for the standard model 
because of the smaller proportion of massive stars. In con- 
trast, adopting a metallicity 20% of the solar value would 
increase Wu a by about 20% because metal-poor stars are 
hotter than solar-metallicity stars. We also tested the in- 
fluence of a possible leakage of some ionizing photons from 
H11 regions into the ambient ISM by allowing 50% of the 
line photons produced by stars in the standard model to 
escape without absorption in the Hi envelopes of the birth 
clouds (Oey & Kennicutt 1997, and references therein). 
While the Ha/H/3 ratio was found to be only 12% lower in 
this case than in the standard model, the Ha equivalent 
width increased by 60%. 

It is also interesting to examine the shape of the total 
effective absorption curve defined by equation (4) and its 
dependence on the age of the stellar population. In Fig- 
ure 5, we show r x (t) at different ages for our standard 
model (computed from eqs. [26] and [27]). Initially, for 
t < £bc> i s a P ower law of index n = 0.7. At later 
times, the effective absorption curve steepens markedly at 
wavelengths A < 1200 A (log A < 3.1). It is more similar 
to the initial power law in the wavelength range where /3 
is determined (1200 < A < 2500 A), but steeper again at 
wavelengths A > 4000 A (log A > 3.6). The reason for this 
is that short-lived, massive stars, which produce most of 
the luminosity at A < 1200 A, are embedded in the birth 
clouds, where their radiation is heavily absorbed. Their 
contribution is substantial and nearly constant across the 
wavelength range 1200 < A < 2500 A but less important 
at longer wavelengths. The steepening of f A at A > 4000 A 
is also enhanced by the appearance in the ambient ISM of 
long-lived stars with a strong 4000 A break in their spec- 
tra, which contribute more light redward than blucward 
of the break. The absorption in the birth clouds, there- 
fore, has a major influence on the shape of the effective 
absorption curve. 



It is worth pausing here to emphasize the important 
implications of the finite lifetimes of the birth clouds. As 
mentioned in §1, emission lines appear to be more atten- 
uated than the continuum radiation in starburst galaxies. 
This led Calzetti (1997) to suggest that ultraviolet-bright 
stars must be somehow physically detached from the gas 
they ionize. In our model, the lines produced in Hn regions 
and the non-ionizing continuum radiation from young stars 
are attenuated in the same way by dust in the outer Hi en- 
velopes of the birth clouds and the ambient ISM. However, 
since the birth clouds have finite lifetimes, the non- ionizing 
continuum radiation from stars that live longer than the 
birth clouds is attenuated only by the ambient ISM. The 
ultraviolet and optical continuum radiation, therefore, ap- 
pears to be less attenuated than the lines. This can be seen 
clearly in Figures 3 c and 4 c from the rise in the Ha/H/3 
ratio accompanying the drop in the Ha equivalent width 
at nearly constant ultraviolet spectral slope as the effec- 
tive absorption optical depth of the birth clouds increases 
above unity. Thus, the finite lifetime of the birth clouds 
is a key ingredient to resolve the apparent discrepancy be- 
tween the attenuation of line and continuum photons in 
starburst galaxies. 

Since each parameter in our model has a specific influ- 
ence on the various integrated spectral properties of star- 
burst galaxies, we can gain insight into the physical origin 
of the relations defined by these observed quantities. An 
inspection of Figures 2-4 leads to the following conclu- 
sions. (1) The wide observed range of ultraviolet spectral 
slopes can be most naturally understood as a spread in the 
effective absorption optical depth in the ambient ISM, al- 
though variations in the lifetimes of the birth clouds could 
also play some role. It is worth mentioning that the steep 
ultraviolet spectra of the galaxies with largest [3 cannot 
originate from old stars because these galaxies also have 
large W-a a . (2) As noted above, absorption in the birth 
clouds is required to explain the large Ha/H/3 ratios ob- 
served in galaxies with intermediate ultraviolet spectral 
slopes. (3) The fact that L dust /ii 6 oo, Lh«/£h/3, and (3 
are all very low at one end of the observed relations and 
all very high at the other end further suggests that varia- 
tions in the effective optical depths of the birth clouds and 
the ambient ISM may be related (Figs. 2 and 3). To ex- 
plore the effects of varying simultaneously fy C and fy M , 
we assume that they are related by fy C — 2ry M and fix 
all other parameters at the standard values given by equa- 
tion (30). The results are shown in Figure 6. Evidently, 
the models reproduce well the observed relations between 
£dust/£i600, Lna/Lnp, and (3 in Figures 6a and 66. They 
also fall roughly in the middle of the observed ranges of 
Wiia and (3 in Figure 6 c. 

An important result of our analysis is that the shape 
of the relation between Ldust/^1600 and (3 actually re- 
quires that the exponent of the effective absorption curve 
^ism x ^-n DC n ~ 0.7. In Figure 6, we also show mod- 
els with f^ c = 2f^ SM as above, but for n = 0.4 and 1.0 
instead of 0.7. These models miss the observed relation, 
the sequence being too steep for n — 0.4 and too shallow 
for n = 1.0. This observational constraint on n applies 
essentially to f A SM , since idust/^1600 and (3 are almost in- 
dependent of t®° for Ty° > 1. For n = 0.7, therefore, 
the relations defined by Ldust/£i600, Ln a /Ln0, and [3 can 
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Fig. 5. — Effective absorption curve (as defined by eq. [4]) of our standard model at different effective starburst ages, as indicated. 
The dashed curve is from Calzetti et al. (1994). For reference, the wavelength range where the ultraviolet slope /3 is determined 
and the wavelengths of the Ha and H/3 lines are indicated at the bottom. All curves are normalized to unity at 5500 A. 
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be understood most simply as a sequence in the overall 



dust content of galaxies. Our assumption of f- 



BC 



2f. 



ISM 



should be taken merely as illustrative. With n = 0.7, the 
observations can be matched by a variety of models that 



range from f 



BC 



-ISM 



= at one end of the relations 
1.5 at the other 



ISM 



(with (3 = -2.3) to f^ c > 1.0 and f, 
end (with (3 w +0.2). For reference, the fraction of the 
total stellar radiation absorbed by dust in the models of 
Figure 6 ranges from 0% to about 90% along this sequence 
and is roughly 70% for our standard model. Non-ionizing 
continuum photons account for about 90% of the heating 
in all models except for those with very low dust content, 
in which the relative contribution by La photons is more 
important (over 20% for (3 < —2.2). As Figures 2-4 in- 
dicate, the intrinsic scatter about the relations at fixed (3 
can arise from variations in the effective age of the star- 
burst, the parameters controlling the birth clouds, and the 
fraction of dust in the Hn regions (and possibly the IMF 
and metallicity of the stars). 

Our models arc also consistent with observations of 
near-infrared H-recombination lines in starburst galaxies. 
The fluxes of the Paschen (3 (1.28 /im) and Brackett 7 
(2.17 /im) lines of 11 galaxies in the Meurer et al. (1999) 
sample are available from Calzetti et al. (1996). We com- 
pute the P/3 and B7 luminosities, Lpp and Lb 7 , using 
the same procedure as described in §2 for Lhq and -Lh/3- 
Figure 7 shows the H/3/B7 and P/3/B7 ratios as a func- 
tion of the Ha/H/3 ratio for the observations and for the 
same models as in Figure 6. Sequences of models are 
straight lines in these log-log diagrams, reflecting the as- 
sumed power-law form of the effective absorption curves 
in the birth clouds and the ambient ISM. Figure 7 shows 
that n = 0.7 is consistent with, although not tightly con- 
strained by, the observed strengths of optical and near- 
infrared recombinations lines in nearby starburst galaxies. 
As noted above, the exponent n of the effective absorption 
curve is much better constrained by the relation between 
the ratio of far-infrared to ultraviolet luminosities and the 
ultraviolet spectral slope. 

It is worth noting that the observed mean relations can 
be reproduced by the following simple recipe: use an effec- 
tive absorption curve proportional to A~ 0,7 to attenuate 
the non-ionizing continuum radiation and emission lines 
(as predicted by dust-free case B recombination) from each 
stellar generation, and lower the normalization of the curve 
typically by a factor of 3 after 10 7 yr to account for the dis- 
persal of the birth clouds. This recipe approximates the 
results of our full calculations shown in Figures 6 and 7 but 
differs from them in that it omits the absorption of ionizing 
photons by dust in H11 regions (i.e., it assumes / = 0.0). 
As mentioned earlier, this has no effect on the ultraviolet 
spectrum and a negligible effect on the far-infrared lumi- 
nosity (< 1% underestimate) and the Ha/H/3 ratio (from 
< 1% to about 7% overestimate along the observed rela- 
tions). The effect is more important for the Ha equivalent 
width (from 3% to over 70% overestimate). Of course, the 
simple recipe above does not allow one to reproduce the 
scatter in the observed relations, which requires variations 



in several parameters of our model. However, if only the 
mean properties are considered, the simple recipe should 
suffice for many purposes. Both our complete model and 
the simple recipe derived from it can be incorporated eas- 
ily into any population synthesis model. 

In an attempt to understand the physical basis for the 
observed relation between idust/^1600 and /?, we plot other 
properties of the galaxies against (3 in Figure 8: the oxygen 
abundance of H11 regions O/H and the absolute B mag- 
nitude Mb- The oxygen abundances of 31 galaxies of the 
Meurer et al. (1999) sample, measured in the same aper- 
ture as (3, are available from Calzetti et al. (1994), while 
the absolute B magnitudes of all the galaxies are available 
from Kinney et al. (1993). The Spearman rank correlation 
coefficients are r s = 0.76 for the relation between O/H and 
(3 and r s = —0.50 for that between Mb and (3. For the dif- 
ferent numbers of galaxies in each case, these correlations 
are significant at the 4<j and 3tr levels. Thus, they are 
probably both statistically significant, especially the one 
between O/H and (3 (which was noted by Calzetti et al. 
1994 from a slightly different sample), but both have very 
large scatter. This suggests that neither O/H nor Mb are 
entirely responsible for the relations between -kdust/£i600> 
Lyia/L}if3, and {3. The ultimate physical explanation for 
these relations therefore requires further investigation. 

For reference, we have also computed models using the 
Calzetti et al. (1994) effective absorption curve. This was 
derived from the analysis of IUE spectra of a sample of 
33 starburst galaxies, nearly all of which are included in 
the Meurer et al. (1999) sample. Calzetti et al. (1994) 
condensed all the spectra in their sample into six "tem- 
plates" corresponding to different values of the Ha/H/3 
ratio. They then derived an effective absorption curve 
by dividing the five spectral templates with the highest 
Ljia/LHp by the template with the lowest Ln a /Lnp (as- 
suming that the latter represented the dust-free case) . The 
result is strictly an effective absorption curve in the sense 
defined by equation (4), although it is often referred to 
by other names (see footnote 2). Figure 5 shows that the 
Calzetti et al. (1994) effective absorption curve is greyer 
than those derived here. In Figure 6, we show the re- 
sults obtained when adopting the Calzetti et al. (1994) 
curve for fj^ M and assuming fy C = for consistency. 7 
As anticipated, the models predict a steeper relation be- 
tween Ldust/^1600 and (3 than is observed. The greyness 
of the Calzetti et al. (1994) curve probably results from 
a bias inherent in the method used to derive it. As Fig- 
ure 66 shows, galaxies with the lowest observed Ha/H/3 
ratios extend over a fairly wide range of ultraviolet spec- 
tral slopes. Thus, the template spectrum with the lowest 
Laa/Ln/3, built by Calzetti et al. (1994) from eight spec- 
tra with —2.08 < (3 < —1.33, may be relatively but not 
perfectly dust free, as is assumed by the method. 

4. COMPARISONS WITH OBSERVATIONS: THE SPATIAL 
DISTRIBUTION OF DUST 

Up to now, we have compared our model with obser- 
vations only in terms of the effective absorption curve, 
i.e., without specifying the spatial distribution and op- 



7 To extend the Calzetti et al. (1994) curve into the near-infrared over the range 0.8 < A < 2.2 (im, we use the formula provided by Calzetti 
(1999). To compute r{ SM at wavelengths 912 < A < 1200 A, we extrapolate the curve linearly in log A. This has a negligible influence on 
the results because typically less than 25% of Ldust is produced by the absorption blueward of 1200 A. In fact, adopting f^ SM = fjfoo f° r 
A < 1200 A gives similar results. 
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Fig. 7. — (a) H/3/B7 ratio and (b) P/3/B7 ratio plotted against Ha/H/3 ratio. The data points are from Calzetti et al. (1996) 
and include 11 galaxies from the Meurer et al. (1999) sample (filled squares) plus two additional nearby starburst galaxies (open 
circles). The solid and dotted curves are the same models as in Fig. 6. Typical measurement errors are indicated in each panel. 



tical properties of the dust. In this section, we study 
the influence of different spatial distributions of dust on 
the integrated spectral properties of galaxies. As shown 
in §3, the effective absorption curve favored by observa- 
tions of nearby starburst galaxies is significantly shallower 
than most known extinction curves. This is a natural 
consequence of any spatial distribution of dust in which 
the least extincted stars contribute most to the emergent 
light at short wavelengths. Such an effect has been quan- 
tified for many distributions of dust, both patchy and 
smooth (e.g., Natta & Panagia 1984; Caplan & Deharveng 
1986; Bruzual, Magris, & Calvet 1988; Disney , Davies, & 
Phillipps 1989; Witt et al. 1992; Calzetti et al. 1994, 
1996; Puxley & Brand 1994; Di Bartolomeo, Barbara, & 
Pcrinotto 1995; Bianchi, Ferrara, & Giovanardi 1996; Witt 
& Gordon 1996; Gordon et al. 1997). These studies, how- 
ever, did not include the finite lifetimes of stellar birth 
clouds. Thus, they cannot account for the different atten- 
uation of line and continuum photons in starburst galaxies. 
We now investigate the influence of the spatial distribution 
of the dust on the effective absorption curve in the context 
of our model. 

For simplicity, we assume that the optical depth of the 
dust has a power-law form 

tx cx \- m , (31) 
where the choice of m is discussed below. This is consistent 
with recent evidence that the optical depth of the dust in 
nearby starburst galaxies does not exhibit a strong feature 
near 2175 A (Gordon et al. 1997). We further assume 
for simplicity that rf c has the same wavelength depen- 



dence as t\. For comparison, Figure 9 shows t\ plotted 
against A for Milky Way-type, LMC (30 Doradus)-type, 
and SMC-type dust and for the two cases of forward-only 
and isotropic scattering (eq. [19]). These optical proper- 
ties are based on the standard Draine & Lee (1984) grain 
model but with the proportions of graphite and silicates 
adjusted so as to fit the observed mean extinction curves 
of the three galaxies (Pei 1992). The inclusion of scatter- 
ing makes t\ greyer because the albedo is larger at optical 
than at ultraviolet wavelengths, the effect increasing from 
Milky Way- to LMC- to SMC-type dust. Power-law fits 
to these curves over the wavelength range 0.1 < A < 1 /im 
yield values of m between 1.1 and 1.5, with m = 1.3 rep- 
resenting a middle value (Fig. 9). 

The ambient ISM of starburst galaxies cannot be ide- 
alized as a foreground screen if the dust within them is 
similar to that in the Milky Way, the LMC, or the SMC 
because the effective absorption curve would then take the 
form fp M cx A~ m with 1.1 < m < 1.5. As shown in §3, 
this would lead to a relation between Ldust/^1600 an d (i 
shallower than observed. We now examine the other two 
models of the ambient ISM specified in §2. Again, these 
should be regarded only as convenient approximations to 
more realistic distributions of dust. In the following, we 
adopt the same effective starburst age, t, lifetime of birth 
clouds, tBG, and fraction of dust in the ionized gas, /, as 
in our standard model in §3 (eq. [30]). 

We first consider the mixed slab model for the ambi- 
ent ISM. In this case, f^ SM is related to the optical depth 
through the slab in the direction normal to the surface, 
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Fig. 8. — (a) Oxygen abundance plotted against ultraviolet spectral slope for galaxies in the Meurer et al. (1999) sample discussed 
in §3. The abundance measurements are from Calzetti et al. (1994). For reference, the solar abundance is (0/H) Q » 8.9. (b) 
Absolute B magnitude plotted against ultraviolet spectral slope for the galaxies in the Meurer et al. (1999) sample. The adopted 
Hubble constant is Ho = 70 kms _1 Mpc _1 . Typical measurement errors are indicated in each panel. 




Fig. 9. — Optical depth of graphite-silicates dust plotted against wavelength. The dotted lines are the true absorption curves 
of the Milky Way, LMC, and SMC, while the dashed lines are the corresponding curves when isotropic scattering is included 
(eq. [19]). These are based on the Draine & Lee (1984) model but with the proportions of graphite and silicates adjusted so as 
to fit the observed mean extinction curves of the three galaxies (Pei 1992). The solid line is a power law of index m = 1.3. All 
curves are normalized to unity at 5500 A. 
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Fig. 10. — Effective absorption curve (as defined by eq. [4]) of a model with effective starburst age t — 3 x 10 s yr in which the 
ambient ISM is represented by a mixed slab (solid line). The optical depth of the dust is assumed to have the form t\ cx A~ m 
with m = 1.3 (dotted line). The short-dashed curve is the standard model of eq. (30), and the long-dashed curve is from Calzetti 
et al. (1994). All curves are normalized to unity at 5500 A. 



cx A~ m , by equations (10) and (23). By analogy with 
our standard model in §3, we adopt fy C = 1.0 for the opti- 
cal depth of the birth clouds and take Ty — 0.5. Figure 10 
shows the total effective absorption curve at t = 3 x 10 8 yr 
with to = 1.3. A comparison with Figure 5 shows that 
this is very similar to the total effective absorption curve 
of the standard model with n = 0.7 in §3. To understand 
the origin of this result, we plotted f^ SM as a function of 
t^ 1 and noticed that the simple formula f^ SM « (j^) 1 ^ 2 ap- 
proximates fj; with an accuracy or 5% or better over the 
whole range 1 < < 10 (at low optical depths, f^ SM is 
closer to rjj 1 ). This weak dependence of f^ SM on rjj 1 arises 
from the fact that, for large rjj 1 , sources located at high 
optical depths from the surface of the slab contribute very 
little to the emergent light. In the ultraviolet, therefore, 
where r|' is greater than unity, the wavelength dependence 
of f* SM is closer to \~ m l 2 than to A _m in the model of 
Figure 10. Hence, for reasonable assumptions about the 
optical properties of the dust, a mixed slab model for the 
ambient ISM can produce an effective absorption curve 
similar to that required by the observed relation between 
£dust/£i600 and (3. 

Mixed slab models, however, cannot accommodate the 
full range of observations. In Figure 11, we show the ana- 
log of Figure 6 for models with Ty C = 2iy and m = 1.3. 
At low optical depths, i.e. for fy C < 1, the models re- 
produce well the observed relation between -Ldust/-ki600 
and j3, as noted above. At larger optical depths, however, 
the models deviate markedly from the observations. This 
upturn, which persists over the range 1.0 < to < 1.5, is 
entirely a consequence of absorption in the ambient ISM 
because, for large optical depths of the birth clouds, i.e. 
for Ty C > 1, (3 and £dust/£i600 are controlled by f^ SM 
(Figs. 2c and 2d). The reason for this behavior can be 



traced back to the fact that fj[ SM scales only as (rf) 1 ^ 2 for 
> 1. As Ty increases in Figure 11, approaches a 

power law of index to/2 over a range that expands toward 
progressively longer wavelengths (since cx TyX~ m ). 
Once the power-law shape extends over the entire ultravio- 
let wavelength range, increasing Ty no longer affects (3 but 
continues to alter l/dust/^ieooj Lyi a /Lnp, and Wn a . This 
is an intrinsic limitation of the mixed slab model. For large 
Ty, (3 cannot reach the very large values observed because 
sources at high optical depths from the surface of the slab 
contribute very little to the emergent ultraviolet radiation. 

We now consider the discrete cloud model for the am- 
bient ISM. In this case, f^ SM is related to the optical 
depth per cloud, cx A~ m , and the mean number of 
clouds along a line of sight, n, by the expression f^ SM = 
n [1 — cxp(— tJ)] (eqs. [10] and [25]). Thus, n does not af- 
fect the wavelength dependence of f^ SM . To illustrate the 
effective absorption curves of the discrete cloud models, we 
fix Ty C = 1.0 and experiment with different optical depths 
of the clouds, t v — 0.1, 0.3, and 0.5. By analogy with our 
standard model in §3, we adjust n so that the mean opti- 
cal depth along a line of sight is tTt v = 0.5 in each case. 
Figure 12 shows the total effective absorption curves at 
t = 3 x 10 s yr for these models with m = 1.3. The mod- 
els exhibit the characteristic signatures of the absorption 
in the birth clouds (Fig. 5). As t v increases, f A becomes 
greyer because the wavelength below which the clouds ab- 
sorb all incident radiation, i.e. where fj~ M « n, becomes 
progressively larger. Figure 12 indicates that the model 
with t v = 0.3 has the ultraviolet slope closest to that 
of the standard model in §3. For reasonable assumptions 
about the optical properties of the dust, therefore, dis- 
crete cloud models also lead to effective absorption curves 
similar to that required by the observed relation between 
idust/ii600 and (3. 
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Fig. 11. — Sequence of models in which the ambient ISM is represented by a mixed slab with different dust content and f v 



BC 



The optical depth of the dust is assumed to have the form t\ oc \~ m with m = 1.3. Crosses mark the models with fy ~ 0.1, 0.5, 
1.0, 2.0, and 5.0 (in order of increasing /?) to indicate the scale. The data points in (a), (6), and (c) are the same as in Figs. 2, 3, 
and 4, respectively. 




Fig. 12. — Effective absorption curves (as defined by eq. [4]) of models with effective starburst age t — 3 x 10 8 yr in which 
the ambient ISM is represented by a Poisson distribution of clouds. The optical depth of the dust is assumed to have the form 
t\ oc A~ m with m = 1.3 (dotted line). The models are shown for different cloud optical depths, Ty — 0.1, 0.3, and 0.5 (as 
indicated). The short-dashed curve is the standard model of eq. (30), and the long-dashed curve is from Calzetti et al. (1994). 
All curves are normalized to unity at 5500 A. 
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Fig. 13. — Sequences of models in which the ambient ISM is represented by a Poisson distribution of clouds with different dust 
content and fy C = 2nr v . The optical depth of the dust is assumed to have the form t\ oc A~ m with m = 1.3. The models are 
shown for different cloud optical depths, t v = 0.1, 0.3, and 0.5 (as indicated). Dotted lines join models with Ty C = 0.5, 1.0, 
and 2.0 (in order of increasing 0) to indicate the scale. The data points in (a), (b), and (c) are the same as in Figs. 2, 3, and 4, 
respectively. 
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The advantage of discrete cloud models over mixed slab 
models is that they can account for the full range of ob- 
servations. In Figure 13, we show the analog of Figure 6 
for models with Ty C = 2nT v , n a continuous variable, 
Ty = 0.1, 0.3, and 0.5, and m = 1.3. The models with 
t v = 0.3 agree with all the observations. In particu- 
lar, they match remarkably well the relations between 
idust/^1600, Ln a /Lft{3, and (3 over the entire range ob- 
served. Models with t v = 0.5 and 0.1 lie above and below 
those with t v = 0.3, as expected from the different effec- 
tive absorption curves in Figure 12. The mean number of 
clouds in the models with t v = 0.3 ranges from n = 
at one end of the relations (with (3 = —2.3) ton « 5 
at the other end (with (3 rs +0.2). These values of n are 
small enough that the distribution of dust can be regarded 
as patchy over most of the range. We also find that the 
relations between £dust/£i600> £h q /£h/3, and (3 are very 
stable to m over the range 1.0 < m < 1.5 for t v = 0.3, 
but curiously less stable for t v — 0.1 and 0.5. The above 
values of t v would be somewhat smaller if we had adopted 
an optical depth t\ with a strong feature near 2175 A. 

Our results demonstrate that, for reasonable assump- 
tions about the wavelength dependence of the dust optical 
depth t\ , discrete cloud models can account for all the ob- 
served integrated spectral properties of nearby starburst 
galaxies. Poisson distributions of clouds have traditionally 
been used to interpret variations in the observed redden- 
ing of stars at the same distance within the Milky Way. 
These observations can be well represented by a model 
with two types of clouds, one with Ay k> 0.2 and another, 
less abundant one with A v w 0.9 (Spitzer 1978). This is 
similar to our model, which has a population of optically 
thin clouds in the ambient ISM and a population of opti- 
cally thick birth clouds. It is worth emphasizing that in 
any model with discrete clouds, the optical depth has a 
large variance along different lines of sight within a given 
galaxy, depending on the number of clouds encountered. 
In particular, some lines of sight can be quite opaque, even 
for T v = 0.3 and n ~ 1. For starburst galaxies with the 
largest values of £dust/£i600) n ~ 5 is indicated, and most 
lines of sight are opaque. Thus, there is no contradiction 
between the relatively small values of t v found here and 
the very dusty appearance of some starburst galaxies. 

5. CONCLUSIONS 

We have developed a simple model to compute the ef- 
fects of dust on the integrated spectral properties of galax- 
ies, based on an idealized prescription of the main fea- 
tures of the ISM. Our model includes the ionization of 
Hii regions in the interiors of the dense clouds in which 
stars form. Emission lines from Hii regions and the non- 
ionizing continuum from young stars are attenuated in the 
same way by dust in the outer Hi envelopes of the birth 
clouds and the ambient ISM. However, since the model also 
includes the finite lifetimes of the birth clouds, the non- 
ionizing continuum radiation from stars that live longer 
than the birth clouds is attenuated only by the ambient 
ISM. We show that this can fully resolve the apparent dis- 
crepancy between the attenuation of line and continuum 
photons in starburst galaxies. This enables us, in turn, to 
interpret in a consistent way all the observations of a ho- 
mogeneous sample of nearby ultraviolet-selected starburst 



galaxies, including the ratio of far-infrared to ultraviolet 
luminosities (I/dust/-^i60o), the ratio of Ha to H/3 lumi- 
nosities (Lhu/Lhp), the Ha equivalent width (Wh q ), and 
the ultraviolet spectral slope (/?). 

The different parameters in our model, including the 
effective age of the starburst, the lifetime and effective op- 
tical depth of the stellar birth clouds, the effective optical 
depth in the ambient ISM, and the fraction of dust in the 
ionized gas, each have a specific influence on the integrated 
spectral properties L dust /L wm , Lua/Lnp, W Ha , and (3. 
This provides new insights into the origin of the mean re- 
lations defined by the data and the scatter about these 
relations. In particular, the relation between the ratio of 
far-infrared to ultraviolet luminosities and the ultraviolet 
spectral slope in starburst galaxies reflects the wavelength 
dependence of the effective absorption in the ambient ISM. 
We find that a power law of the form f^ SM cx A -0 7 ac- 
counts remarkably well for all the observations. The re- 
lation between £dust/£i600 and (3 can then be interpreted 
as a sequence in the overall dust content of the galax- 
ies. Interestingly, this relation is accompanied by much 
weaker trends of the oxygen abundance and the optical 
luminosity with the ultraviolet spectral slope. The fact 
that our model reproduces the observed spectral proper- 
ties of nearby starburst galaxies relatively easily leads us 
to suspect that, with suitable adjustment of the parame- 
ters, it could also reproduce those of more quiescent (but 
still star-forming) galaxies. 

The effective absorption curve required by the obser- 
vations is much greyer than would be produced by a fore- 
ground screen of dust like that in the Milky Way, the LMC, 
or the SMC. We have explored whether this could be ac- 
counted for by a much steeper wavelength dependence of 
the optical depth, i.e. t\ cx A~ m with 1.0 < m < 1.5, 
combined with a more realistic spatial distribution of the 
dust. We find that a mixed slab model for the ambient 
ISM can produce the required effective absorption curve 
for low dust content but cannot explain the observations of 
starburst galaxies with very reddened ultraviolet spectra. 
In contrast, we show that a random distribution of dis- 
crete clouds provides a consistent interpretation of all the 
observed integrated spectral properties of starburst galax- 
ies. While these results were anticipated in some previous 
studies, we have shown here for the first time how to rec- 
oncile them with the large Ha/H/3 ratios and other obser- 
vations. We also find that the optical depths of the clouds 
favored by our analysis are similar to those inferred from 
the statistics of stellar reddening in the Milky Way. 

The model we have developed for computing the absorp- 
tion of starlight by dust in galaxies can be combined easily, 
by design, with any population synthesis model. The ob- 
served mean relations for starburst galaxies can also be re- 
produced by the following simple recipe: use an effective 
absorption curve proportional to A -0 7 to attenuate the 
line and continuum radiation from each stellar generation, 
and lower the normalization of the curve typically by a fac- 
tor of 3 after 10 7 yr to account for the dispersal of the birth 
clouds. This recipe accounts at least as well as the one by 
Calzetti et al. (1994, and as modified by Calzetti 1997, 
1999) for the effects of dust on the non-ionizing contin- 
uum radiation. In addition, it fully resolves the apparent 
discrepancy between the attenuation of line and contin- 
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uum photons in starburst galaxies. We believe, therefore, 
that our model and the recipe derived from it provide sim- 
ple yet versatile tools to interpret the integrated spectral 
properties of starburst and possibly other types of galax- 
ies. In future work, we plan to apply them to the growing 
body of observations of high-rcdshift galaxies. 
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